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Abstract We discuss the general theoretical arguments advanced earlier for the 
T = global phase diagram of antiferromagnetic Kondo lattice systems, distin- 
guishing between the established and the conjectured. In addition to the well- 
known phase of a paramagnetic metal with a "large" Fermi surface (Pl), there 



■ is also an antiferromagnetic phase with a "small" Fermi surface (AFs). We pro- 



vide the details of the derivation of a quantum non-linear sigma-model (QNLaM) 
representation of the Kondo lattice Hamiltonian, which leads to an effective field 
theory containing both low-energy fermions in the vicinity of a Fermi surface and 
, low-energy bosons near zero momentum. An asymptotically exact analysis of this 

' effective field theory is made possible through the development of a renormaliza- 

O ! tion group procedure for mixed fermion-boson systems. Considerations on how 

O . to connect the AFs and Pl phases lead to a global phase diagram, which not only 

puts into perspective the theory of local quantum criticality for antiferromagnetic 
, heavy fermion metals, but also provides the basis to understand the surprising re- 

^ ' cent experiments in chemically-doped as well as pressurized YbRh2Si2. We point 

OO I out that the AFs phase still occurs for the case of an equal number of spin-1 /2 lo- 

^^O ■ cal moments and conduction electrons. This observation raises the prospect for a 

22 ', global phase diagram of heavy fermion systems in the Kondo-insulator regime. Fi- 

. nally, we discuss the connection between the Kondo breakdown physics discussed 
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here for the Kondo lattice systems and the non-Fermi liquid behavior recently 
studied from a holographic perspective. 
PACS numbers: 

1 Introduction 

Critical phenomenon associated with a second-order phase transition is formu- 
lated in terms of Landau's notion of spontaneous symmetry breaking. Consider 
an antiferromagnetic transition of a rotationally-invariant spin system. Above the 
Neel transition temperature, T > T^, spins retain the rotational invariance. Below 
the Neel temperature, T < T^, the rotational invariance is spontaneously broken, 
leading to an antiferromagnetic (AF) state in which the order parameter - the stag- 
gered magnetization, m - acquires a finite value. Approaching the critical point, 
the order parameter vanishes but strong spatial fluctuations of the order parameter 
remain. Indeed, the criticality is characterized by this fluctuating order parameter, 
m(x), and is described by a Ginzburg-Landau theory in c/-spatial dimensions. 

What happens when the Neel temperature is tuned to zero? Hertz formulated 
the problem in the context of itinerant systems, in which the order parameter m 
is associated with a microscopic magnetic state known as the spin density wave 
(SDW)™. Hertz further assumed that the fluctuations of this SDW order parameter 
are the only critical mode. For any quantum statistical system with Hamiltonian//, 
its partition function, Z = Tre~'^^, where j3 = 1 /kT, can be considered as a quan- 
tum mechanical evolution operator along imaginary time — ;t, where T is defined 
in the range (OjfzjS) with an associated periodic boundary condition. When 7)v is 
nonzero, finite-size effects eventually cut off the fluctuations along the imaginary 
time axis leaving classical critical behavior associated with the spatial fluctuations 
of the order parameter. When 7}v — ?> 0, on the other hand, the imaginary time axis 
has an infinite extent. Correspondingly, the critical theory is a quantum Ginzburg- 
Landau theory of m(x, t) in d +z dimensions, where z is the dynamic exponent. 
The temporal fluctuations keep track of the quantum nature of the collective fluc- 
tuations, but the degree of freedom itself remains to be m, the coarse-grained 
staggered magnetization, which is a classical variable. 

In recent years, it has become possible to test this picture of order-parameter 
fluctuations with the explicit observation of the magnetic quantum critical points 
(QCPs) in heavy fermion metals^i^^SiJ^. At the same time, a large body of theo- 
retical work has emerged on the quantum phase transitions in the Kondo lattice 
Hamiltonian. 

The purpose of this paper is multi-fold. We will first discuss the early theoret- 
ical and experimental motivations for the local quantum criticality with a Kondo 
breakdown. We describe the aspects in which this type of QCP differs from the 
Landau picture of order-parameter fluctuations, as well as the extensive theoretical 
developments and the pertinent experimental results. The emphasis here is how the 
well-established heavy fermion state, a paramagnetic metal with a "large" Fermi 
surface (Pl), can be critically broken down. 

We proceed to consider the antiferromagnetically ordered part of the Kondo- 
lattice phase diagram. Our recent work— introduced a quantum non-linear sigma- 
model (QNLcjM) representation of the Kondo lattice Hamiltonian. An asymptot- 
ically exact renormalization-group (RG) analysis establishes the existence of an 



3 



antiferromagnetic phase with a "small" Fermi surface (AFs). Here, we provide 
the details of the derivation of this QNLcjM representation, and describe the RG 
procedure for this mixed fermion-boson system. 

These considerations have lead to a proposal for a global phase diagra m 
This phase diagram puts the studies of local quantum criticality into a larger per- 
spective. In addition, recent experimental work on the Ir- and Co- doped YbRh2Si2— 
have uncovered some surprising varieties of quantum phase transitions in heavy 
fermion metals. Related observations have appeared in pure YbRh2Si2 under pres- 
sure^-- and Ge-doped YbRh2Si2 These recent theoretical and experimental works 
have opened up a new direction for the study of quantum criticality and novel 
phases in heavy fermion metals. 

We also consider the case of heavy fermion systems in the regime of Kondo in- 
sulators. We show that the small-Fermi-surface antiferromagnetic phase, AFs, also 
exists in the l-l-l filling case. Correspondingly, we propose that tuning a Kondo in- 
sulator towards larger RKKY interactions can induce it into an antiferromagnetic 
or paramagnetic metal with a small Fermi surface. 

The focus of this article will be on the theoretical issues. More discussions on 
the experimental systems, with more references, can be found in a recent article^ 
as well as in some of the comprehensive review articles^ii^Sii^. 



2 Local quantum criticality and Kondo breakdown from the paramagnetic 
side 

2.1 Kondo breakdown and local quantum criticality 

Early indications for the failure of the Hertz picture came from inelastic neutron 
scattering measurements of the dynamical spin susceptibility ;if (q, ft)). Following 
initial studies of dynamical scaling in heavy fermion metals^, measurements of 
% (q, ft)) in the heavy fermion compound CeCu6_xAux at the critical concentration 
jfc w 0. 1, show that the order-parameter fluctuations contain an anomalous critical 
exponent, and display the property of ft)/r scalin g^^'"*^'^^ . Both are properties of 
an interacting fixed point. By contrast, the SDW QCP descriptio n ^^i^'i-^^ of such 
systems would have an effective dimensionality of + z > 4, the upper critical 
dimension of the quantum Ginzburg-Landau theory, and would correspond to a 
Gaussian fixed point with essentially mean-field exponents and violation of (o/T 
scaling. 

These results motivated the theoretical proposal of a new class of QCPs. The 
key characteristic is the emergence of critical modes which are inherently quantum 
mechanical. For the Kondo lattice system at hand, these modes are associated with 
a critical Kondo breakdown. This Kondo breakdown picture was already being 
considered prior to these experimental developments in some related theoretical 
models-^^i^2ii2. in light of the observed interacting behavior, speculations were put 
forward— that the Kondo breakdown underlies the interacting behavior observed 
in the dynamical spin susceptibility. Soon thereafter, concrete theoretical formula- 
tions were advanced^^*i2j^. The Kondo breakdown picture has subsequently been 
discussed in related formulation s '^^i^^ . 
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To put things in perspective, the traditional view of heavy-fermion metals is 
that the ground state is a Kondo singlet. This is a singlet state formed among 
all the local moments and all the conduction electrons. This Kondo singlet state 
does not involve any spontaneous symmetry breaking, but involves a macroscopic 
order. This is a quantum entangled state. A continuous breakdown of the Kondo 
entanglement gives rise to the quantum critical modes. 

A series of work s 54i55^23^63^72,64^58^2i ,73,22 jj^yg studied the Kondo breakdown 

effect within an extended dynamical mean- field theory (EDMFT)^^i^ of the 
Kondo lattice Hamiltonian. The EDMFT equations have been studied using e- 
expansion RG procedure and various numerical methods. When the Kondo break- 
down occurs at the magnetic QCP, the criticality depends on not only the order- 
parameter fluctuations but also the new Kondo-breakdown critical modes. Corre- 
spondingly, fractional critical exponents and (o/T scaling arise, in a way that is 
consistent with the experimentally observed dynamical spin susceptibility. 

This local quantum critical picture also has implications for electronic excita- 
tions. Kondo resonances carry both spin (1/2) and charge (e). A Kondo breakdown 
at the magnetic QCP leads to a jump of the Fermi surface at the magnetic quan- 
tum critical poin t^'^i^i^'^i'^i^^ , as well as a Kondo-breakdown energy scale, 
that continuously goes to zero at the magnetic QCP. Both properties also operate 
when the Kondo breakdown occurs away from the magnetic Q(;;;p,57j54,55j49j38 

The Fermi surface jump has been probed extensively in several heavy-fermion 
metals. In YbRh2Si2, it has been demonstrated that Hall coefficient is dominated 
by the normal component and therefore probes the Fermi surface^. At low tem- 
peratures, the Hall coefficient displays a rapid crossover as a function of the non- 
thermal control parameter (a relatively small magnetic field). This crossover ex- 
trapolates to a jump at the magnetic QCP in the limit of zero temperature, pro- 
viding evidence for a jump of the Fermi surface. In CeRhlnj^^i^, evidence for a 
Fermi surface jump accompanied by a mass divergence has come from measure- 
ments of the de Haas-van Alphen effect.^. 

The evidence for a Kondo-breakdown energy scale going to zero at the mag- 
netic quantum critical point has come from the measurements of both the Hall 
effect— and thermodynamic quantities^. 



2.2 Extended dynamical mean field theory of Kondo lattice 

The collapse of Kondo effect has been extensively studied from the paramag- 
netic side. A number of effects come into play: the RKKY interactions promote 
magnetic order, the Kondo interactions favor Kondo-singlet formation, and the dy- 
namical competition between these two types of interactions are important for the 
transition. One suitable microscopic approach is the extended dynamical mean- 
field theory (EDMFT) ^^i^^i'^ . Here we outline the basic equations and the two type 
of solutions 54,55,23,72,63,2i,73.22 _ ^ niore detailed summary can be found in Ref.— . 
We consider the Kondo lattice model: 

Here, = Xlkcj EkV'kcj Vkcj describes a band of free conduction c— electrons, with 
a bandwidth VK. M'k = Lc^a:S, -Scj specifies the Kondo interaction of strength 
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Jk', here the conduction electron spin Sc,; = ^Lcjct' i'^oa'Wa' where T is the 
vector of Pauli matrices. Finally, ,3ff = j ^/yS, • Sj is the magnetic Hamiltonian 

for the spin-5 /—moments, S,. The strength of the exchange interactions, /,,, is 
characterized by, say, the nearest neighbor value, /. 

Within EDMFT, the dynamical spin susceptibility and the conduction-electron 
Green's function respectively have the forms xiH: = [^q+M{(0)]~ \ and G(k, e) = 
[e + ji — e]i — E{e)]^^ . The correlation functions, ;t(q, ft>) and G(k, e), are momentum- 
dependent. At the same time, the irreducible quantities, M{(0) and Z(e), are momentum- 
independent. They are determined in terms of a Bose-Fermi Kondo model, 



^£S- i^p + ^lp) +Y^wp^l^p . (2) 



p 



The fermionic (Cpa) and bosonic (^p) baths are determined by self-consistency 
conditions, which manifest the translational invariance, Xioc{(0) = HqXi^.j^'^)' 
and Gioc{£) = LkG(k,£). The (0+ 1) -dimensional quantum impurity problem, 
Eq. (|2]l, has the following Dyson equations: M{o)) = Xo^{Oi) + l/Zioc(w) and 
Z{e) = Go He) - l/Gioc(e), where Xo'M = -g^Lp2wp/{co^ -wD and Go(e) = 
l/(£ — Ep) are the Weiss fields. The EDMFT formulation allows us to study 
different degrees of quantum fluctuations as manifested in the spatial dimension- 
ality of these fluctuations. The case of two-dimensional magnetic fluctuations are 
represented in terms of the RKKY density of states that has a non-zero value at 
the lower edge, e.g., 

p/(x) =£5(A--/q) = (1/27)0(7- |x|) , (3) 
q 

where is the Heaviside step function. Likewise, three-dimensional magnetic 
fluctuations are described in terms of a pi {x) which vanishes at the lower edge in 
a square -root fashion, e.g.. 



Piix) = {2/nf)Vl^-x^0{I-\x\). (4) 

The reduction of the Kondo-singlet amplitude by the dynamical effects of 
RKKY interactions among the local moments has been considered in some de- 
tail in a number of studies based on EDMF T 54i55,23,72,63,2i ^73^22 _ irrespective of 
the spatial dimensionality, this weakening of the Kondo effect is seen through the 
reduction of the E^^ scale. 

Two classes of solutions emerge depending on whether this Kondo breakdown 
scale vanishes at the AF QCR In the case of Eq. (|4]l, TJj* ^ has not yet been com- 
pletely suppressed to zero when the AF QCR dc, is reached from the paramagnetic 
sideQThe quantum critical behavior, at energies below TJj*^, falls within the SDW 
jypg25j3ij33 -pjjg zero-temperature dynamical spin susceptibility has the following 
form: 

Xifl,(o) (5) 
/(q) - laco 

' However, it can still go to zero inside the AF region— i—i^. 
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Here /(q) =lq Jq, and is generically (q — Q)^ as the wavevector q approaches 
the AF ordering wavevector Q. The QCP is described by a Gaussian fixed point. 
At non-zero temperatures, a dangerously irrelevant operator invalidates the (o/T 
scaling-^ii^. 

Another class of solutions corresponds to E^^^ = already at dc- It arises in 
the case of Eq. O]), where the quantum critical magnetic fluctuations are strong 
enough to suppress the Kondo effect. The solution to the local spin susceptibility 
has the form 

^^''^''^ = m+A{-Lrw{co/T) - 

This expression was derived within EDMFT studies, through the aid of an £- 
expansion approach to the Bose-Fermi Kondo model. At the AF QCP, the Kondo 
effect itself is critically destroyed. The calculation of the critical exponent a is 
beyond the reach of the e-expansion. In the Ising-anisotropic case, numerical cal- 
culations have found a w 9 7 23^21^73^22 

The breakdown of the Kondo effect not only affects magnetic dynamics, but 
also influences the single-electron excitations. As the QCP is approached from the 
paramagnetic side, the quasi-particle residue zl °= {t>*)'^, where b* is the strength 
of a pole in the conduction-electron self-energy L, goes to zero. The large Fermi 
surface turns critical. 



3 Quantum non-linear sigma model representation of the Kondo lattice 

The local quantum criticality discussed in the previous section implies a break- 
down of the Kondo effect at the magnetic QCP. In order to explore this issue 
further, we study the Kondo effect deep inside the antiferromagnetic region of 
the phase diagram. This regime arises when the Kondo coupling is infinitesimal 
compared to the AF exchange interaction among the local moments, i.e. in the 
regime 

7a- < / < W (7) 

where W is the width of the conduction electron band. 

In the regime specified by Eq. (|7]l, our strategy is to expand the Kondo lattice 
with respect to the limit of Jk = 0, where the local moments and conduction elec- 
trons decouple. The appropriate approach to implement this strategy comprises 
two steps. First, we introduce a representation of the Kondo lattice Hamiltonian 
in terms of a Quantum Nonlinear Sigma Model (QNLaM). This representation, 
which results in an effective low-energy field theory containing mixed fermions 
and bosons, has been introduced in Ref.^; the goal of this section is to provide 
the details of the derivation. Second, this effective low-energy theory is analyzed 
asymptotically exactly in terms of an RG method; we will go through this in some 
detail in the following sections. Related studies have also been pursued in Ref.— . 
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3.1 Summary of the mapping 

Since an explicit demonstration of the mapping will take a fair amount of space, 
in this subsection we summarize the essential points. Subsequent subsections will 
provide the details. 

The Kondo lattice Hamiltonian has been specified in Eq. ([l]). For now, we will 
consider the electron concentration, x per site, to be such that the Fermi surface 
of J^fc alone does not touch the antiferromagnetic zone boundary. Later, we will 
discuss the modifications necessary for the more general case. 

We expand around the limit = 0, where the local-moment and conduction- 
electron components are decoupled. We will consider, for simplicity, square or 
cubic lattices, although our results will be generally valid provided the ground 
state is a collinear antiferromagnet. J^f can be mapped to a QNLcjM by standard 
means-Sli^. The low-lying excitations are concentrated in momentum space near 
q = Q (the staggered magnetization) and near q = (the total magnetization being 
conserved): 

2S/ ^ T]xn(x, t)^1 - {2adL{x, t))^ + 2fl^L(x, t) (8) 

where x labels the position, T]x = ±1 on even and odd sites, a is the lattice constant, 
and we have used S = 1/2. The linear coupling n ■ Sc cannot connect two points 
on the Fermi surface and is hence unimportant for low energy physics (such a 
kinematic constraint has appeared in other contexts, e.g. Ref.— ); see Fig.[T^. The 
Kondo coupling is then replaced by an effective one, S-Sc~^ fl^'L ■ s^, correspond- 
ing to forward scattering for the conduction electrons; see Fig.[T^. 




Fig. 1 (a) When the Fermi surface of the conduction electrons does not intersect the boundary 
of the antiferromagnetic Brillouin zone (AFBZ, the dashed lines), the only component of the 
Kondo coupling to the QNLctM fields involves forward scattering for the conduction electrons 
near the Fermi surface, (b) (from Ref.--) When the Fermi surface of the conduction electrons 
does not intersect the AFBZ, the Kondo coupling connects the re-diagonalized fermions (whose 
Fermi surfaces are given by the solid lines) to the QNLcfM fields. This coupling is also forward 
scattering in the AFBZ. 
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The mapping to the QNLaM can now be implemented by integrating out the 
L field. The effective action is 



QNLcjM ' 



QNLcjM 



— / dPxdx 

2g 



^erry 



(9) 



='k I d^'xdx [Sf (x, t) ■ (p(x, t) 



■^c = j d''Kd£Y,¥l{^,i£){i£~^K)V<y{1^,i£)+^^ I 

where = Vf{K — Kp). The Berry phase term for the n field, oS^Berry, is not im- 
portant inside the Neel phase, which is very different from the ferromagnetic case. 
What is meant by "the Berry phase" requires some clarification. Certainly some 
aspects of the geometric term do indeed contribute to the physics, but this will 
be spelled out in the next section. We have introduced a vector boson field (p 
which is shorthand for n x The n field satisfies the constraint = 1, which 

is solved by n = (tt, a), where n labels the Goldstone magnons and a = — 
is the massive field. We will consider the case of a spherical Fermi surface; since 
only forward scattering is important, our results will apply for more complicated 
Fermi-surface geometries. The parameters for the QNLaM will be considered as 
phenomenological- , though they can be explicitly written in terms of the micro- 
scopic parameters. The effective Kondo coupling is /I = Uj^/ {Adia''), which will 
be explicitly demonstrated below. 

This summarizes the structure and setup of the effective field theory for the 
antiferromagnetic phase of the Kondo lattice model. We now describe the details 
of how this is done. 



3.2 Coherent state representation of the partition function 

We set up our notation by considering in some detail the standard case of the 
Heisenberg model. This will help us perform the analogous mapping for the Kondo 
lattice model, which is essentially the same but includes conduction electron cou- 
pling. We will focus on a square lattice with nearest-neighbor (nn) and next- 
nearest-neighbor (nnn) spin-exchange interactions, I\ and I2 respectively. 

nn nnn 

H = ^lijSi ■ Sy = /i ^ S,- ■ Sj + h S, ■ Sj = /i ^ Sx ■ Sjt+a + /2 ^ Sx ■ Sj-+j3 
'■j {ij) (((/)) 

(10) 

where x runs over all lattice sites, a runs over the nn sites for each lattice site x, 
and j3 runs over nnn sites for each site x. After the mapping is completed, it will 
be clear that the coupling constant g of the QNLaM can be tuned by changing /] 
and 12- 

The coherent state spin path integral representation of the partition function for 
quantum spin systems is now a standard formalism that can be found in textbooks 
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(e.g.^^^). We will briefly sketch the main idea. The partition function can be 
written 



„ M-l 

Z = Tre-P"=lim h^A(T) Yl{A(Tk+i)\l - eH\A(T,)) (11) 



where the many-particle basis is a direct product of single-site spin coherent states: 
|A(t)) = nvl'^.v(f))- Here M is the number of discrete time slices in the Trotter 
decomposition, and £ = P/M. Recall that at each lattice site, (i2|S|i2) = S£2, 
where the unit spin vector is represented by 

Q = =S w 8/5 = (cos sin 0, sin sin 0, cos 0). 

y'S(S + 1) 

Using this (overcomplete) basis, the matrix elements of the Heisenberg Hamilto- 
nian can thus be written, to leading order in e, 

{A{Zk+i)\H\A{Zk)) = S^{A{Tt+i)\A{Tt)) 

(12) 

This allows us to write the Hamiltonian in terms of classical variables Qx- To 
linear order in e, 

(A(T,+i)|(l-e//)|A(T,)) = {A{Zk+,)\A{Zk))e-'"^'^'^^ (13) 

where 

H,,{t) = S^hY.iixiT) ■ iix+aiT) +S\'£^xiT) ■ iix+pi^) (W) 

x.a x,p 

The penalty for the classical representation is the additional overlap (A {Xk+i) \A{Tk)) 
which is the Berry phase accumulated from the adiabatic evolution from the time- 
slice Tji to Tk+i- Including the Berry phase accounts for quantum corrections, and 
is crucial to obtain the proper mapping to the QNLaM. We can write it more 
clearly as follows: 

(A(-r^+i)|A(Ti)) = ll{Qx{Tk+i)\£ix'{^k)) 

x,x! 

^W{Qx{Tk+x)\^x{Tk)) 

X 
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In the partition function we need an infinite product of such overlaps, which leads 
to a continuum representation in imaginary time 



M-l 



lim Yl{A{T,+ i)\A{z,))= hmlle 



-/S[l-cose,(T^)][0,-{Tt+l)-0^{Tt)] 



k=l 



lim e 



-''SLv.*[l-COsev(TiO][0A-{Tj,+ l)-0x(Ti)] 



^e-'-SLv.ff«"t[i-cose,(T)];* 
_ ^-/5E,m(I2,) 



(16) 



where (o{Qx) = Jo '^'^[^ ~ '^^^ 0v(t^)]^ is the Berry phase for a single spin at 
site X. Note that the total Berry phase contribution to the action in the path inte- 
gral is given by the sum of the Berry phases of all the lattice site spins: ,5^ = 
iS'£,x('^{^x)- We have represented the Berry phase of a single spin with a set of 
parameters 6^ and (j)^ for familiarity, but we will work with an alternative and 
convenient representation given by 



(o{Q) = / dx I du 

Jo Jo 



du 



dx 



(17) 



where by convention Q{t,u — \) = and Q{t,u = 0) = (0,0,1) = |0) = 
15, m = S) (see, for example, Ref.— ). Cl{x^u) is a homotopically equivalent con- 
tinuous deformation of i2(T). 

Now, in a similar way to what was done above for the Berry phase, we can 
take the continuum limit to express the Hamiltonian term, Hd, as an integration 
over imaginary time. The partition function then becomes 



Z = 



dxHciix) 
dx I du 







du 



dx 



x,p 



(18) 
(19) 

(20) 
(21) 



3.3 QNLaM mapping for the Heisenberg model 

We assume an antiferromagnetic order in which case we can represent each spin 
as the sum of a staggered component representing the local Neel field, and 
uniform {q w 0) fluctuations L j, 



fl,(T) = 77,n,(T)^/l-(yL,(T)l +yL,(T) (22) 



a' 
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The factor 77^ G ±1 depending on which sublattice x falls in, while the factor 
a'' /S 'm front of the uniform fluctuation field L ensures that integrating around any 
small volume will yield the total magnetization contained in that volume. Here a 
is the lattice constant. Recall that the spin variable is constrained by the condition 
Six ■ Six = 1 at each site. With the above choice, this constraint now becomes 
n.v • = 1 and • = 0. Note that the total number of degrees of freedom in 
the system remains the same because in the n, L representation we must restrict 
ourselves to the antiferromagnetic Brillouin zone (AFBZ); there are twice as many 
degrees of freedom on half as many sites. 

We now wish to write the action in terms of n and L rather than CI . We will 
consider the Berry phase first, then the He] term. 



3.4 Berry phase 

Let us first consider what happens when we substitute this expression for the spin 
into the Berry phase part of the action. We need the expressions for the T and u 
derivatives: 



-rr— = i2„ = 777n„ L-L„ n + — L„ 

au 76 S 

—-=£l, = nYn,-^ — {L-U)n+—U 

OT 76 S 



(23) 
(24) 



where we temporarily dropped the site index, x, and instead use subscripts to de- 



note differentiation. We have also defined 7 = y 1 — ] ■ Plugging this into 
equation ( l20l ): 



JO 



777n+^L ) • ( 7727V, xn^- 77^ ^(L-L„)n„xn 



+T7yn« xLT-772 — (L-L„)nxnT + ^^^(L-L„)(L-LT)nxn 

77 a^"^ a'^ 77 a^'^ cP^ 

(L • L„)n X Lt + 7] 7— L„ X (L ■ Lt)L„ x n + -^-L,, x L 

76^ S Y ^ ^ 



To simplify this equation requires knowing that n„, and L are all perpendicular 
to n. This means their triple product must vanish: L ■ n„ x = 0. We also neglect 
terms higher than linear order in L [terms quadratic in L are small compared to 
those kept in Eq. (l34l)1. leading to 

7] Vn • n« X + 77 V-^rn • ("« x + L„ x 



/■/3 /-l 

Yb » iSY / dT du 

xJo Jq 



(26) 
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Note also that r]^ = I and 77^ = rj^. We then obtain 



du 



« iSY f dz [ 
"Jo Jo 

= iSY dx f 
^ Jo Jo 

= iSY dx du 
V-/o Jo 



du 



T]n • n„ X lit + — n • (n„ x L-;- + L„ x n^-) 



ad d f dn, 
S du \ dx 



dux driTc a'-' 
T?.vnt • ^ X — + — n, 
du dx S 



dWx ^ '^L.v ^ 
du dx du dx 



dUx a!^ d ( dw^x 

T].,n, • — X — + — — 11;, • — X L, 
du dx S dx \ du 



(27) 



In the second line we have restored the full notation, while the third line can be 
written with total derivatives since the terms proportional to cancel thanks 
to the triple product identity 1^ • b x c = — b • 1^ x c. The second term in the third 
line vanishes after integrating the total x derivative and using the periodicity of 
the fields. The third term in the third line can be integrated over u, and the value 
at M = is zero due to the orthogonality at the north pole. We finally find, 



yg = iSYrjx dx du 
^ Jo Jo 



dux dux 

du dx 



r Jo 



C?T Lv • llv X 



dnx 
dx 



(28) 



The first term is precisely the Berry phase for the Neel component n, while the 
second term is something additional that must be added to the total action. Al- 
though both terms came from the expression for the Berry phase of fl, it is only 
the first term that is often referred to as the Berry phase for the antiferromagnet. 
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3.5 Hamiltonian 



Next we compute the contribution to the action from the HamiUonian Hci{t) ex- 
pressed in terms of n and L fields. Plugging ( l22l ) into ([21]), 



77.vn,(T) W 1 - — L;,(t) + — L,(t) 



T]^+„n,+a(T)^/ 1 - ( yL,+ „(T) I + yL,+ „(T) 



2 d 

a" 



T],,n,(T)^/ 1 - ( yL,(T)') + yL^(T) 



'7;.+/3n.v+/3(T)\/ 1 - -^Lv+i3(T) + -^L;.+/3(t) 



2 d 



/ 1 fa'' 
J7.n.(T) l--f-L,(T; 



-L,(t) 



1 fa! 



2 \ S 



Lv+0!('C) 



1 /fl' 



-5^/2 1 



■L,(t) 



-Lv+/3(t) 



(29) 



Since n is a unit vector, we have the identity n;r • n,, = 1 — ^ (n.v — n,.) • (iiji — n,.). We 
also know that at every site iiv ■ Lv = 0. Using these two identities and dropping 
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the T label for brevity, the Hamiltonian can be expressed as follows: 



l-2K-"v)-(n.v-nj,) 



^2d 



1 :^i^x + ^x+a) 



a'^ a'' 



a'' a'' 

2d 



2 , T 2 



+ 



a 



252 



(30) 



where we used 2La • Lv+a = L? +L2^„ — (L^ — Lv+a)^ and similarly for /3. We 
have written the expression in this way in order to take advantage of a Taylor 

series expansion between different lattice sites: — n^K, a{a ■ V)«^ ^ and 

n^^^ — Rj fl\/2(i3 • H . The index b runs over the components of the 

vector field. Note that nearest neighbor (nn) sites x and x + a are separated by a 
distance a, while next nearest neighbor (nnn) sites x and x + j3 are separated by a 
distance a\/2. Expressing all lattice differences in this way leads to: 



Hci = S^h £ { r]xTlx+c 



x,a 

^\ 
252 . 



l--[(a-V)«;f 



fcl2 



,2d 



+- 



a 

"252 L 

a^+y2 
52 



( 1 - rixTlx+p ) (Li + L2+^ ) - 2«2 [()3 . v)L^] 



il2 



(31) 
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After summing over lattices sites and neighbors, each term proportional to a''^^ 
sums to zero, so. 



Hcl = S^h £ <^ T];tT?:t+o 



2d 



l-y[(«-V)«?]2 



a 
"252 



(1 - t],t],+„)(l2 + L?+ J ~ a^[{a ■ V)L^f 



^2d 



a 

252 



(32) 



Now, to our order of approximation, L^^^ w L^. Also, since a runs over nearest- 
neighbors, while j8 runs over next-nearest-neighbors, we have i^x'Hx+a = — 1 and 
'7aT?.v+/3 = +1- On the square lattice, the number of nn and nnn sites is 2d, and to 
avoid double counting we divide by 2. 



252 



2d' 

2 



4L2-^(VL,)2 



-S^dl^ll-'MWn,) 



2a^ 



(VL,)2 



(33) 



Note that (Vn^)^ should be interpreted as L^^i LfeLi ■ We also used L^.a [(a ■ 

V)n*]2 = + ^ + = I.f(Vn,)2, and similarly for j3 and L terms. 

To be consistent with our expansion in small powers of a and 1 /S we should 
also ignore gradient terms like (VL)2. The expression simplifies to 



H,i =S^dY^\h-h + ^(/i - 2/2)(Vn, 



S^d^,,eil2- h) +S^^ih -2l2)'£0^n,f +2dha^'£^l (34) 



Finally, we take the continuum limit with the correspondence E ^ — s- a / d x, 

Hci{t) = consti+S^a^~''{Ii/2~l2) J d''x{yn{x,T)f +2dlia'' j d''xL^{x,T) 

(35) 

We have introduced the constant factor const i = jVsitedS^{l2 — h ) which is unim- 
portant for our purposes. 
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3.6 Completing the square and the QNLaM mapping 
At this point the total action is given by 



Z = I S)n^L 5(n2 - l)5(n • L)^-'^!"^^! (36) 
y[n,L] = comfi+^^[n]+5V-^(/i/2-/2)^'^<iT j d'' x {Vn{x, T)f 
+2dha'' dT I d''xV-{x,T) 

-i f^^dx I d'^x (^L(x,t) ■n(x,T) x ^^j^^ (37) 



^i[n] EE /5£77, r dT f du U ■ ^ X 

r Jo Jo V OU dT 



(38) 



where the Berry phase o5^g[n] now only depends on the n field. The delta func- 
tionals enforce the local constraints. To deal with them, we use the integral repre- 
sentation of the delta functional and introduce a scalar Lagrange multiplier field 
a(x,T): 



(39) 



It is now clear that the functional integral is Gaussian with respect to L, so we 
may "complete the square" and integrate it out completely. If we use the identity 
/^Le/-?L'+bL ^ a:'^/2(det^)-i/2gbV4f ^ tjie correspondence is ^ = 2dha'^ and 
b = m X n — ian. We also need the quadruple vector product identity (n x h)^ = 
j,2jj2 _ . jj^2^ ^j^j jj^g relations h • n = and = I and n • n x h = 0. This leads 



to 



4? 



(h + a ), and hence 



y[n,a] 



(40) 



y[n,a] = /3comri - -log— — ^+5V-^(/i/2-/2) / dx / d'x{Vn{x,T)f 
I 2Ai\a" Jo J 



1 



Sdliad Jo 



dT / d"x 



dn{x,z) 
dx 



const2 + [n] + ^ t [ d' 
2g Jo J 

I rP r 
+— / dx / d''xa^(x,x) 
IgcJo J ^ ' 



+ a^(x,T) 



(Vn(x,T)r + 



1 /<9n(x,T)^ ^ 



dx 



(41) 
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where we have defined 



c = 2aSh Vd^^-^-j^ (42) 
1 

const2 = j3 consti ~ ^ log ^ ^ (44) 
2 2fl/ifl" 



Recall that M is the number of time slices and may be considered to be of order 

All that remains is to perform the gaussian integral over a: / i^a e ' = 
^lgc%^ = e{i/2)iog(2.?cjr'^)^ Quj. fljial answer becomes: 



.2 l)e-nn] 



Z= l^nd{n- 

c fP r 
y[n] = constj + ^i[n] + — / dT d'^x 

2g Jo J 



, 1 / dn(x, t) 



(45) 



(46) 



where constj, = constj - jlog(2gc7z'^) = P^KiteS^dih - h) - log(27r''^). These 
results for the constants c and g agree with— who considered first, second and 
third neighbor couplings. Clearly, by adjusting Ii and /j we can tune g. This model 
can also be expressed in terms of the spin-wave stiffness and transverse magnetic 
susceptibility. They are given by 



= -=5V-''(/i-2/2) 

1 1 



eg Aa'^lid 



(47) 
(48) 



Notice that Xi- is independent of h, which means that nnn interactions do not 
renormalize the transverse magnetic susceptibility. 



3.7 QNLaM mapping for the Kondo lattice model 



We have now shown how to map the quantum Heisenberg AF to the QNLaM. 
Next, we want to incorporate the Kondo interaction which couples the local mo- 
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ment spin, S, to the conduction election spin, Sc- This adds the following term; 

JKS'£iix{T)-ScA^)=JKSa-'' j d''xdTSc{x,T)-n{x,T) 

T].,(n(x,T) •Sc(x,t)) 



JKSa'"^ / d'^xdT 



fad \^ a'^ 

'^y^^xY^^^'^^J +yL(x,T)-s,(x,T) 

kJkJ t/''xt/TL(x,T) •s,(x,t) (49) 

The last line follows due to n ■ w 0. The latter is because we have, as discussed 
earlier, chosen to work with a Fermi surface that does not intersect the AFBZ 
boundary (see Fig.[T^); what remains of the Kondo interaction is the (nearly) for- 
ward scattering channel for the conduction electrons. The assumption we make 
is not necessarily that the density of conduction electrons is infinitesimally small, 
but only that it does not intersect the AFBZ boundary. Specifically, we require 
Q > IKp- In section 132] we will discuss the modifications to the theory when 
the Fermi surface does indeed intersect the magnetic zone boundary—, i.e. when 
Q < 2Kf. 

Let us return to the action for the quantum AF before completing the square 
(equation 137b. and add to that the above Kondo coupling. The total action for the 
Kondo Lattice Model, J^klm^ now has something extra coupled to the L field: 



Z= y^n^L^a^v^"^^v^5(n2-l) 

X e-'5"A-LM[n,L,a,s,]-.r, [y/t , 1//]-.^^ [n] ^-^q^ 

^jCLM[n,L,a,Sc] = consti + S^a^^'' (h /2 - h) dx j d'' x {Vn{x, T)f 
+2dh a!^ 1^ dT j d'^xl? (x, t) 

- y dx j d''xL{x, t) • ( m(x, t) x — JkSc(x, t) 

-ja(x, T)n(x, t) I (51) 



where const\ and ^g[n] are as defined previously, and S^c[w\w\ is the con- 
duction electron component of the action. Just like before, the functional inte- 
gral is Gaussian with respect to the L field. With the identity / ^Le '^"^^ = 
^ / det i^e'' , the correspondence is now C, = Idlia'' and b = m xh — JKSc — 
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ian. The important quantity is: 



4C Sdha" 



2„2 



— (n X h) — UkSc - nxh+an-nxh — UkSc -0x11 + 7^s 



2„2 



'UkOCSc • n + an • n X h + i'/^-Sc • n - a n 



(52) 



To simplify this expression, we need the following identities: n^ = 1, h • n = 0, 
(n X h)^ = h^n^ — (h • n)^ = — h^, n • n x h = 0, Sc • n w 0, and 



a.a.p.y.S \ / 



This leads to: 



4C 8t//ifl'' 



37; 



- - 2i7xSc •nxh+— ^ £v/^v/^- 2vA,,VtV^|v^; - aV 



Note that the terms that came from serve only to renormalize the direct quadratic 
and quartic fermion couplings, which can be incorporated into J/'dY^, i//]. So after 
integrating out the L field we find: 



Z = f ^n^a^v^^^v^5(n2-l)e-^^«["'«'*-l-^^[v'''V]-^i[n] 



(Vn(x,T))2 + ^' ^ ^ 



+A j d''xdT ^Sc(x, t) • n(x, t) x ^ 
: dx 1 d''x a2(x,T) 



5t 



_1_ /-/^ 
2gc. 



(53) 



The constants c, and const2 are defined exactly as before, and the new Kondo 
coupling constant is: 



~ 4t//ia'' 



(54) 
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Finally, we integrate out the a field which only contributes the same constant as 
before. The final result is: 



^QNLcjM[n] = const-i 



2g 7o 



1 / ^n(x,T) 

c2 V dx 



(Vn(x,T))2 



^A:[n,Sc] = A j d''xdT ^Sc(x, t) •n(x, t) x — 
^c[V^"^</] = / ^/''^^/e£4(K,/e)(/e-^j,)v/a(K,ie) + a J / (58) 



(55) 



(56) 



(57) 



Note that terms from have been absorbed in ii and ^k- This completes the map- 
ping from the microscopic Kondo Lattice Hamiltonian to the effective field theory, 
as claimed earlier. 

Now that we have demonstrated this mapping in detail, we are confronted 
with executing the RG analysis. The next section is devoted to this theoretical 
development. 



4 Renormalization group analysis and the antiferromagnetic phase with 
Kondo breakdown 

In the previous section we presented the details on the construction of the effective 
field theory corresponding to the antiferromagnetic Kondo (Heisenberg) lattice. 
The end result is a QNLaM coupled to itinerant electrons with a Fermi surface. 
An RG analysis of this field theory is much more involved than field theories 
typically encountered in high energy physics. These complications stem from the 
existence of a Fermi surface^^. We wish to treat the gapless fermions and bosons 
on an equal footing instead of integrating out the fermions first; doing the latter 
would introduce undesirable non-analyticities, as first pointed out byi^^ii. 

We have developed an RG approach for mixed fermion-boson theories^, build- 
ing on the method for purely-fermionic problems The approach is not techni- 
cally difficult, but the subtleties of the scaling procedure are not straightforward 
either. Our method meshes the scaling of bosons, which takes place along all di- 
rections in the momentum space with respect to a point in momentum space, and 
that of fermions, which involves only one direction locally perpendicular to the 
Fermi surface. This is illustrated in Fig. |2l For a long exposition of this scaling 
procedure, we refer the reader to our recent paper-Si. Here, we wish to emphasize 
only one particularly pertinent aspect. Our mixed fermion-boson approach is sim- 
plified in the antiferromagnetic case considered here by the fact that the bosonic 
modes of the QNLaM action [Eq. ( I56H has a dynamic exponent z = 1 . 
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2Ab 




Ab « 
Af « Kf 



Fig. 2 Kinematics of the RG rescalingi22 for a system containing botli gapless bosons, whose 
momenta are centered around a point in the momentum space, and gapless fermions, whose 
momenta are confined to the vicinity of a Fermi surface. Ab and Af are the cut-off momenta for 
the bosons and fermions, respectively, a is the lattice constant in real space. 



4.1 The case of Fermi surface not intersecting the antiferromagnetic zone 



The analysis based on the combined fermion-boson RG has been discussed else- 
where^l. The effective Kondo coupling, expressed in the QNLaM representation 
by Eq. ( 157b . is marginal at both the tree, one-loop, and infinite-loop levels. The 
tree level and one-loop analyses directly use the procedure of Ref.™, while the 
infinite-loop analysis was based on a decomposition of the Fermi surface into 
patches. 

Here we present an alternative infinite-loop analysis without appealing to patch- 
ing arguments. The goal of this analysis is to establish a type of Migdal's Theorem 
which states that the tree-level result is the entire story. Since we found marginality 
at the tree-level, this is the exact answer to all orders in the limit where A /Kf — )■ 0. 

For concreteness, we consider a spherical Fermi surface although generaliza- 
tions to non-nested generic Fermi surfaces can be readily made. For a spherical 
Fermi surface we can write the Fermi momentum integral in terms of spherical 
coordinates: 



boundary 




(59) 



The most relevant part of the above is 




(60) 
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Now the kinetic part of the fermions can be written, 

=^F~^ J dkd'^-^QKde - Vf fc) (61) 
We define new dimensionless variables: 

e=Ae (62) 

k = Ak (63) 

Qk = (64) 

/s:^-'aV'> = V^^'V (65) 

Note that the angular components of fermionic momenta are untouched. We now 
have: 



= Jdkd''-'^QKdexif\ie-VF~k)x}f 



(66) 



The important difference from the aforementioned patching argument is that now 
the fermionic fields contain factors of Kf. Plugging this into the Kondo coupling 
we find (note that the QNLaM rescaling is identical to what was done in Ref.— ) 

= /s:;-^Ai+i+''+i-3+i-(^+3)/2;i,^y dkd'^~^Q.KdMUdS)[\ii'''^f{G)7t)l61) 

rfZi— ^± / dW'^ ^KdM'' qdcbl^f^ -^{m)] (68) 

For c/ = 2 we have 



In the non-spin-flip channel: 



J- _ ^\-dy^l + l+d+\+d+l-3+l 



-2((rf+3)/2) ^_ J ^i^d- 1 QKd£d''qida)id'^q2dan [ij/^ v/((y TTTT )] 



^?L, / dkd'^~''Q.Kd£d'^qid(bid'^q2d6yi[^'^-^{mn)] (70) 

F ■' 



For d = 2, 

Kf Na 

We have thus shown that for any li > 1 the Kondo vertex will have associated 
with it positive powers of 1 /N/^ ■ Because of this, as the number of powers of 
Jk increases, so does the suppression factor 1 /N^ ■ The only exception is for a 
series of diagrams corresponding to a chain of particle-hole bubbles in the spin- 
flip channel-^. Because the poles are located on the same side of the real axis, 
they make no contribution to the beta functio n Thus, the tree-level result is 
the whole story, and the Kondo coupling is exactly marginal. 

The RG result is further corroborated by a large-N calculation for the conduc- 
tion electron Green's function, which does not contain a pole and, correspond- 
ingly, the Fermi surface is small—. 
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4.2 The case of Fermi surface intersecting the antiferromagnetic Brillouin zone 
boundary 

We now turn to the case where the Fermi surface intersects the AFBZ boundary, 
hi this case, the linear couphng n • between the local moments and conduction 
electron spin cannot be neglected. Until now, we have only considered the term 
L ■ Sc because our assumption has been that the Fermi surface does not intersect 
the AFBZ boundary, i.e. Q > 2Kp. See Fig[T^ and the comments following equa- 
tion ( 1491 ). When Q < 2Kf, the conduction electrons see the AF order parameter 
of the local moments as a staggered scattering potential, resulting in a reconstruc- 
tion of their Fermi surface. The hot spots of the Fermi surface therefore become 
gapped out, as shown in Fig. [ij). The effective coupling"^ between the recon- 
structed quasiparticles and the n field involves a coherence factor containing an 
additional factor of q (measured w.r.t. Q). This linear-momentum suppression fac- 
tor survives beyond the mean-field treatment of the conduction electron band, as 
dictated by Adler's Theorem. Indeed, it can be viewed as a kinematic suppression 
similar to the deformation potential problem of the electron-phonon system—. 

From the RG perspective, the additional factor of q represents a decrease in 
the dimension of the coupling which has the same effect as a derivative coupling: 
Se ■ — qa^ aa-Uq. Now the conduction electron spin is coupled directly to 
which has dimension [Uq] = —d — 1, but the additional factor of q brings the di- 
mension to [q Wq] = —d which has the same value as the vector field [(p,^] = —d 
we considered earlier for the case where the Fermi surface does not intersect 
the AFBZ boundary-^. Therefore, our previous result on the marginality of the 
Kondo coupling is not spoiled when the Fermi surface intersects the magnetic 
zone boundary. 

The marginal nature of the Kondo coupling implies that the effective Kondo 
coupling does not flow towards strong coupling. Correspondingly, there is no 
Kondo singlet formation in the ground state, and no Kondo resonances in the 
single-electron excitation spectrum. The absence of the Kondo resonance implies 
that the local moments remain charge-neutral, and the Fermi surface is determined 
by the conduction electrons alone; such a Fermi surface is called small, and the 
antiferromagnetic phase is named AFs. This is to be contrasted with what hap- 
pens in an antiferromagnetic phase in the presence of Kondo resonances, in which 
the Fermi surfaces are specified by the hybridizing conduction electrons and the 
itinerant /-electrons that describe the Kondo resonances; such a Fermi surface is 
called large, and this antiferromagnetic phase is named AFl. Note that the pres- 
ence of antiferromagnetic order can be incorporated in these definitions through 
an appropriate symmetry-breaking field — . 



5 Global phase diagram of the antiferromagnetic Kondo lattice 

The establishment of the AFs phase, along with the heavy-fermion Pl phase-^*^, 
provide two anchoring points in the heavy-fermion phase diagram at T = 0. In 
order to connect these two phases, Ref.-^ introduced a two-parameter phase di- 
agram, (7^,G), which is shown in Fig. [3^. In units of conduction electron band- 
width, Jk is the Kondo coupling while G corresponds to magnetic frustration or re- 
duced dimensionality. The latter is a measure of the quantum fluctuations amongst 
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the local-moment degrees of freedom. It was also discussed in Ref.— that, for 
sufficiently large G, the conventional AF state of the local-moment component 
will yield to paramagnetic states which either preserve spin-rotational invariance 
(gapped or gapless spin liquid) or break it (spin Peierls). Inspired by the recent 
developments in doped YbRh2Si2, this part of the phase diagram is explicitly in- 
cluded in the global phase diagram~ as shown in Fig. [3]3. Related considerations 
on the global phase diagram are also being made in Ref.—. 



G 




Fig. 3 (a) Global phase diagram with the illustrative Fermi surfaces. (From~.) (b) Global phase 
diagram showing the three different types of trajectories for quantum phase transitions. This 
phase diagram has a remarkable correspondence with what happens in heavy fermion metals. 
(FromiS.) 

This global phase diagram can be described in terms of three possible trajec- 
tories that connect the AFs phase and Pl phase, as specified in Fig.[3j3. Trajectory 
/ goes directly between the AFs and Pl phases, giving rise to a local quantum 
critical point: the collapse of Kondo resonances occurs at the magnetic quantum 
critical point - in the notation of Refs.^'^5^5^, 5c and 5[g^ are located at the same 
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place. As the Kondo effect is continuously broken down at the AF QCP, there 
is a sudden jump between the large Fermi surface and the small one 54,10,55^57 ^ 
and the Kondo-breakdown scale vanishes at the QCP. Correspondingly, the 
residues associated with both the small and large Fermi surfaces vanish as the 
QCP is approached from either sid c^^'^"'^^'^^'^^ . Our considerations of the global 
phase diagram puts local quantum criticality in a larger perspective. 

Along Trajectory //, the transition between the AFs and Pl phases involve an 
intermediate AFl phase. This is the SDW state of the heavy quasiparticles of the 
Pl phase. The magnetic-to-paramagnetic transition is of the SDW typ c^^'-^^'^^ . A 
Kondo breakdown transition can still take place at the AFl — AFs boundary— i^; 
in typical cases, this corresponds to a a Lifshitz transition with a change of Fermi 
surface topology. 

Finally, along Trajectory ///, the transition goes through the intermediate Ps 
phase. This is a paramagnetic phase with a small Fermi surface. The small-to-large 
Fermi surface transition could be either a spin-liquid^2i^ to heavy-Fermi-liquid 
QCP, or a spin-Peierls to heavy-Fermi-liquid QCP~. 

There is overwhelming evidence that pure YbRh2Si2 at ambient pressure dis- 
plays a type-I quantum phase transition. Recently, Friedemann et ah— showed 
that enough Co-doping (of nominally 3% or more), which introduces positive 
chemical pressure, turns the transition into one whose properties are largely com- 
patible with the type-II transition. Preliminary studies have shown that a similar 
effect arises in pure YbRh2Si2 under a sufficiently large pressure^. A relatively 
small amount of (nominally 2.5%) Ir-doping, which introduces negative chemical 
pressure, retains the type-I transition—. By contrast, a larger negative chemical 
pressure, corresponding to a larger amount (nominally 6% or more) Ir- doping, 
turns the transition into one that is compatible with a type-Ill transition. A similar 
behavior has also been seen in Ge-doped YbRh2Si2 ' - . Experimentally, there is 
an indication that the regime corresponding to Ps has non-Fermi liquid behavior, 
raising the exciting possibility that this is in fact a non-Fermi liquid phase. 

6 Kondo insulators 

Heavy fermion metals involve localized magnetic moments and a partially-filled 
conduction-electron band. This can be generically modeled in terms of a Kondo 
lattice Hamiltonian comprising two specifies of electrons: a lattice of spin- 1/2 
local moment, with one per unit cell; and a band of conduction electrons with a 
filling of < Jf < 1 electrons per unit cell. 

When X = 1 , we have an even number of electrons per unit cell and an insulator 
becomes a natural possibility. In fact, the analogue of the Pl phase of the < 
X < \ case is the Kondo insulator—. For x = I, the Kondo-singlet formation in 
the ground state induces delocalized /-electron quasiparticles, which hybridizes 
with the conduction electron band and induces a hybridization gap at the Fermi 
energy. At zero-field, it is traditionally believed that the Kondo insulator is the 
only possible phase. (By contrast, a large magnetic field can suppress the Kondo 
effect altogether, inducing a magnetically ordered metal—.) 

Here, we make the observation that the AFs phase extensively described above 
remains a stable phase for the 1 + 1 "Kondo-insulator filling." The exchange in- 
teraction between the local moments in this regime is generically expected to be 
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antiferromagnetic. In addition, in the parameter regime specified by Eq. (Q, our 
asymptotically exact RG analysis of the Kondo lattice Hamiltonian in its QNLaM 
representation continues to apply, regardless of whether the underlying conduction 
electron Fermi surface intersects the AFBZ boundary. To the same degree of con- 
fidence as in the case of generic filling, the small-Fermi-surface paramagnetic (Ps) 
phase is also expected to occur in the case of 1 + 1 filling when the the parameter G 
(the quantum fluctuations of the local-moment component) becomes sufficiently 
large and for small Kondo coupling. 

The above assumes that there is no perfect nesting of the conduction electron 
Fermi surface. A perfect nesting with x = 1 would only occur for a square lattice 
with strictly no hopping beyond nearest neighbors, and this is unlikely to be the 
case for rare-earth intermetallics. 

These considerations lead to a global phase diagram for Kondo insulators, 
shown in Fig. |4l The list of materials believed to be Kondo insulators is by now 
relatively large. It will be instructive to search for quantum phase transitions out 
of a Kondo insulator and into either the AFs or Ps metallic phases. While it is in 
general hard to predict precisely which trajectory a particular tuning parameter - 
such as pressure or chemical doping - would correspond to in such a global phase 
diagram, what is minimally required is to increase the ratio of the AF RKKY 
interaction to the Kondo coupling. For Ce-based systems, this happens with the 
application of a negative (chemical) pressure. For Yb-based systems, on the other 
hand, this is achieved by applying a positive (chemical) pressure. 




Fig. 4 Global phase diagram for Kondo-insulator systems. The Kondo-insulator phase corre- 
sponds to a Kondo-singlet ground state, and is paramagnetic. AFs describes an antiferromag- 
netic metal state with a small Fermi surface and Ps a paramagnetic metal phase also with a small 
Fermi surface; each has the same meaning as its counterpart in Fig. [3] 
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7 Connection to the holographic theory of non-Fermi liquid 

There has been considerable recent interest in applying the techniques developed 
in string theory to many-body systems. These developments rely on a duality be- 
tween a £/-dimensional strongly-coupled field theory and a quantum-gravity the- 
ory in a weakly curved {d + 1) dimensional anti-de Sitter (AdS^z+i) spacetime. 
Non-Fermi liquid behavior has been studied through a charged black hole in the 
AdS spacetim e 28,29^ 1 1 . 1 5^40 ^j^j^ approach, the fermionic self energy shows a 
marginal-Fermi-liquid-like form—. This has been clarified as originating from 
an infrared fixed point characterized by a (0 + 1) -dimensional conformal field 
theory—, which is dual to a near horizon geometry AdS2 x M''^'. The (0+ 1)- 
dimensional conformal field theory corresponds to a quantum impurity model. 

The factorization of temporal and spatial correlations become more transpar- 
ent in a "semi-holographic" formulation in which a conduction-electron band 
is introduced to be hybridized with the fermions in a strongly coupled field theory 
dual to the gravity theory. In turn, this allows an Anderson-lattice-model inter- 
pretation (see also Ref.— ), in which the fermions of the strongly coupled field 
theory is the analogue of the /-electron degree of freedom of the Anderson lat- 
tice model. The relevance or irrelevance of the hybridization may then be linked 
to the Kondo screening or Kondo breakdown of a class of quantum impurity 
models - the Bose-Fermi Kondo model - coupled to a fermionic bathi^^ii^jZlZi. 
In the Kondo language, a breakdown of the Kondo effect, occurring when the 
hybridization-induced Kondo coupling is irrelevant or marginal in the RG sense, 
leads to a power-law form of the conduction-electron self-energy whose exponent 
is positive; such a form of the self-energy appears in the non-Fermi liquid state of 
the holographic models^^i^^Jiii^iiS. A Kondo-screened state, occurring when the 
hybridization-induced Kondo coupling is relevant in the RG sense, gives rise to a 
pole in the conduction-electron self-energy. No such a form of the self-energy has 
been identified so far in the holographic models. 

These considerations point to a linkage between the Kondo-breakdown physics 
of the Kondo/ Anderson lattice systems and the fate of the non-Fermi liquid state 
dual to the gravity theory with a near-horizon AdS2 x M''" ' geometry. There could 
be at least two possibilities for the latter. One possibility is that the non-Fermi liq- 
uid state is unstable towards an antiferromagnetically-ordered state with a small 
Fermi surface discussed earlier, but retains a similar form in the quantum-critical 
regime; this is the analogue of the Kondo-breakdown physics in the local quantum 
criticality formulatio n ^^i^^ . Note that the entropy vanishes at the local quantum 
critical pointi^,. An alternative is that this state is unstable towards a paramagnetic 
phase with a small Fermi surface^^, as occurring in one formulation of the Kondo 
breakdown-t' J'- . In either case, these considerations suggest that holographic mod- 
els may contain an analogue of the paramagnetic phase with a large Fermi surface, 
whose conduction-electron self-energy singularly depends on energy (containing 
a pole) and smoothly depends on momentum. 

With these considerations in mind, it will be instructive to study itinerant mag- 
netic systems from the gravity side. A recent work has gone along this direction—. 
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8 Conclusions 

Recent theoretical and experimental developments have opened up the issue of 
the global phase diagram in heavy fermion metals. We discussed some of the 
earlier developments that have led to the continued theoretical efforts on the global 
phase diagram, and the striking recent experiments that point to the richness of 
the quantum phase transitions between antiferromagnetic and paramagnetic heavy 
fermion metals. 

We provided the details of the asymptotically exact analysis of the Kondo 
lattice in the antiferromagnetic part of the phase diagram. General considerations 
on the transitions from the Kondo-breakdown antiferromagnetic metal phase with 
a small Fermi surface to the Kondo-screened paramagnetic heavy fermion phase 
have motivated a global phase diagram, which is consistent with earlier studies of 
the Kondo breakdown effect from the paramagnetic side including the distinction 
between the local quantum critical point and the SDW quantum critical point. Our 
considerations also put earlier formulations of the Kondo breakdown effect in the 
paramagnetic region into a general perspective. Model studies that can cover all 
these different phases and transitions in one unified framework are called for. 

We have also discussed the Kondo insulators along a similar line. We have 
proposed a related global phase diagram for such systems, which we hope will 
stimulate future experiments. 

Finally, we have discussed the connection between the Kondo-breakdown lo- 
cal quantum criticality of Kondo lattice systems and the holographic non-Fermi 
liquid behavior. Considerations of magnetic states and quantum magnetic transi- 
tions from a gravity dual are just beginning. Conversely, studies of quantum crit- 
ical behavior in magnetic many-body systems, including heavy-fermion metals 
and insulators, may shed some light on gravity problems. 
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